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Reactions induced by high energy antiprotons on proton on nuclei are accompanied with large 
probability by the emission of a few mesons. Interesting phenomena can be observed and QCD 
tests can be performed, through the detection of one or more mesons. 

The collinear emission from high energy (anti)proton beams of a hard pion or vector meson, can 
be calculated similarly to the emission of a hard photon from an electron [?]. This is a well known 
process in QED, and it is called the "Quasi-Real Electron method", where the incident particle is 
an electron and a hard photon is emitted leaving an ’almost on shell’ electron impinging on the 
target [?]. 

Such process is well known as Initial State Emission (ISR) method of scanning over incident 
energy, and can be used, in the hadron case, to produce different kind of particles in similar 
kinematical conditions. 

In case of emission of a charged light meson, K or p -meson, in proton-proton!anti-proton) colli¬ 
sions, the meson can be deviated in a magnetic field and detected. 

The collinear emission (along the beam direction) of a charged meson may be used to produce 
high energy (anti)neutron beams. This can be very useful to measure the difference of the cross 
sections of (anti)proton and (anti)neutron scattering from the target and may open the way for 
checking sum rules with antiparticles. Hard meson emission allows also to enhance the cross 
section when the energy loss from one of the incident particles lowers the total energy up to the 
mass of a resonance. 

The cross section can be calculated, on the basis of factorized formulas, where the probability 
of emission of the light mesons multiplies the cross section of the sub-process. Multiplicity 
distributions for neutral and charged meson production are also given. 
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Figure 1: Total (black) and elastic (red) cross section for the p + p reaction, as function of antiproton beam 
momentum pp. The contribution of inelastic events is also shown (green line). 


1. Introduction 


Among (anti)proton induced reactions which occur with high probability (the cross section is 
of the order of [mb]), we focus here to light meson production in the kinematical domain accessible 
at the future FAIR facility by the PANDA experiment (the antiproton momentum is in the range 
1.5-15 GeV). 

The world data for the total and elastic cross sections for p+p induced reactions are illustrated 
in Fig. [I] as function of the laboratory antiproton momentum [|T(j|. From the figure we see that 
Gpp/Op'p ~ 1/5 and that the difference between the total and elastic cross section parametrizations 
(green line), which represents the contribution of the inelastic events, is of the order of 40 mb in 
the range accessible at PANDA. 

Figs. 2 shows that the most probable reaction involving pions, corresponds to more than three 
pions in the final state. We parametrize the cross section for the reactions given in Fig. 2 as: Among 
the inelastic processes, the charge exchange process pp —t nh has the largest cross section: 

It is known from QED that two Feynman diagrams contribute to the process of Bhabha- 
scattering: the scattering (t) and the annihilation ( s ) channels. Neglecting the interference one can 
estimate the relative contribution of the scattering diagram I-#/ 2 and the annihilation one ./f /, u | 2 
for the case when the final particles are emitted at large angles in the system of center of mass: 


Il^l 2 _ 1 
L |^| 2 ~ to’ 


(l.i) 


The contribution of scattering mechanism dominates. This statement holds even if hadron form 
factors are taken into account, i.e., one can assume that a similar relative contributions of these two 
mechanisms would take place in case of elastic and inelastic pp —>• pp cross section. The processes 
relevant to "scattering" mechanism are denoted as peripheral processes. 
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Figure 2: Cross sections for different inelastic channels in the p + p reaction as a function of pp\ (left) 
mesonic channels; (right) mesonic and hadronic channels 


In this work we estimate the CE contribution, using the formalism [;8j, [9]. The emission by 
the initial proton of a charged light meson -n or p -meson in proton-proton(anti-proton) collisions 
transforms high energy protons (for example in a proton beam) into neutrons. This effect is ob¬ 
served in accelerator physics [1C]. Note that hard collinear photon emission was also studied in 
Ref. This works uses a frame where the relevant amplitude interferes with the background, 
i.e., all other emission mechanisms, except emission from the electron line, whereas |[i| applies the 
structure function approach (that we use). Although the methods are different, the answers of these 
publications are equivalent for the observables. 

CE reactions may occur in antiproton beams, too. High energy, high intensity antiproton beams 
will be available in next future at PANDA [12], FAIR [|h3|]. Hard n and p meson can be detected 
with high efficiency. Charged meson will be deviated by the 2T magnetic field of the central 
spectrometer, and (anti)neutrons, which are produced at high rate, could be used as a secondary 
high energy beam. We are interested here in ’initial state emission’, as a mechanism to produce a 
(anti)neutron beam with comparable energy of the (anti)proton beam. 

In this work, we propose a description of CE reactions referring to the known QED process 
of emission of a hard real photons by electron (positron) beams at e + e colliders. Such process 
enhances the cross section when the energy loss from one of the incident particles lowers the total 
energy up to the mass of a resonance. This is known as "return to resonance" mechanism. In 
the case of creation of a narrow resonance this mechanism appears through a radiative tail: it is 
the characteristic behavior of the cross section which gradually decreases for energies exceeding 
the resonance mass. This mechanism provides, indeed, an effective method for studying narrow 
resonances like J/'•V. 

For the emission in a narrow cone along the directions of the initial (final) particles, the emis¬ 
sion probability has a logarithmic enhancement, which increases with the energy of the "parent" 
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Figure 3: Feynman diagram for collinear hard photon emission in eT reactions (T stands for any target) 


charged particle. In frame of QED this mechanism is called as "quasi-real electron" mechanism 
(QRE) [§. 

In this work we apply the QRE mechanism to the case of hadrons and, in particular, to the 
collinear emission of a light meson from a (anti)proton beam. We evaluate the cross section for 
this process for single as well as multi pion production, where pions can be neutral or charged. Our 
derivation concerns a kinematical region outside the resonance formation region. 


1.1 Quasi-real electron kinematics with hard photon emission 

Let us consider the radiative process e(p \) + T(p 2 ) —>• e(p\ —k) + y(k) +X (the four momenta 
of the particles arc written in parenthesis), T stay for any nuclear target(proton p, or nucleus A), 
and the final state X is undetected, Fig. [| 

The virtual electron after the hard (collinear) photon emission is almost on mass shell [{8]. This 
property allows to express the matrix element of the radiative process e(p) + T —> e(p — k) + y(k) + 
X in terms of the matrix element of the non-radiative process e + T —>• e + X : 

^yiPhPl) =e^(p 2 ) ^ 1 J i + m £(k)u(pi). ( 1 . 2 ) 

— Zp\K 


In the case when the denominator of the intermediate electi on’s Green function is small \ (p\ — 
k) 2 — m 2 1 -C 2p\pi one can write p\ — k + m = £ v u s (p\ — k)u s (p\ —k) and the matrix element has 
a factorized form. 

The square of the matrix element, summed over the spin states of the photon is: 


II- 


= e 


E l +E U 


m 


CO(E p - CO)(kp) (kp) 2 


II^(P2Mpi-*)f 


(1.3) 


where £ | ■'f{p 2 )u(p\ — k )| 2 is the Born matrix element squared with shifted argument. 

In the case of unpolarized particles, the cross section of the process e{p\) + T (p 2 ) —>• e + 7 + T 
may be written in factorized form: 


dGy(s,x) 

dWy(x) 


do{xs)dWy{x),x = 1 — x, 


a dx 
n x 


[l —x + 


1 


r2n 2 

-x 2 )\n—^--(l-x) 
2 mi 


i x =%; 0 < 00 « 1, — > 1, (1-4) 
t m P 
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Figure 4: Feynman diagram for collinear hard charged pion (p-meson) emission in p(p) + T collisions. 


where E is the energy of the initial electron (center of mass frame implied pi + P 2 = 0), and 
s = {p\+Pi) 2 - 

It is assumed here that the initial electron transforms into an electron with energy fraction 
1 — x and a hard photon with energy fraction x which is emitted within the cone 0 < 0 {) along the 
direction of initial electron. Moreover it is implied that xs > s t h r , where s t / ir is the threshold energy 
of process without photon emission. The logarithmic enhancement originates from the small values 
of the mass of the intermediate electron, which is almost on mass shell. This justifies the name of 
Quasi Real Electron (QRE) method. 

Below we consider a possible extension of the QRE method to the processes with "quasi real" 
(anti)nucleon intermediate state. 

1.2 Application to hadron physics 

Let us apply this formalism to the case of initial high energy proton (anti-proton) beams and 
the emission of a hard pion or vector meson in forward direction, collinear to the beam. We do not 
consider the emission of the pion from the final proton (final state emission), which corresponds to 
a different kinematical region for the intermediate particle. 

For the case of emission of a positive-charged p + or n 1 meson by the high energy (anti)proton, 
the final state consists in a high energy (anti)neutron, accompanied by a positively charged meson. 
The charged meson can be deflected by an external magnetic field, providing the possibility to 
select a high energy neutron beam. In the case of emission of the neutral meson, it can be identified 
by measuring its decay channels. 

Let us consider the reactions (Fig. jdj): 


(1.5) 
( 1 - 6 ) 

5 


p + T —>• n + T + h + 
p + T —>■ n + T + h~ 









Light meson emission 


E. Tomasi-Gustafsson 


where li = p or n and T may be any target (p, n, nucleus..). In QED the interaction is assumed 
to occur through the exchange of a virtual photon. In the present case the particle can be a vector 
meson. However, the nature of the exchanged particle is irrelevant for the present considerations, 
where we focus on small angle charged pion emission and on the factorization properties of the 
cross section. 

The matrix element for collinear Jt(p) emission can be written as: 

^kl{PUP2)) = ^nj{P\ -k : p 2 )^ h P ' l {p]. : P] ~k), 

^l’ T (pi,p 2 )) = JZnT{pi-k,p2)^”{p\,p\-k), ( 1 . 7 ) 

where 

?n & 

^n n {Pl,Pl-k) = 2 — - TUn(Pl~k)Y5U p (pi) 7 

2p\k 

^p n {pi,Pi~k) = 2 8 — -u n (pi-k)eu p (pi). (1.8) 

mj t — 2p\k 

are the matrix elements of the subprocesses: p —>• n + n + and p -X n — p + (or p —>• « + and 

p —»• n + p~). The matrix element of emission of a charged pion is expressed in factorized form 

( |L7| ) of a term connected to the target and a ’universal’ factor for the hard meson emission, which 
generalizes the QED result. Below, we will focus on this last term. 

The relevant cross sections are: 


d (jPT^h+X( s ^ = o nT ^ x {xs)dW h+ {x ), 
da pT ^ h+x (s,x) = o nT ~ rX (xs)dWh_ (x), 
da pT ^ hoX (s,x) = c J pT ^ x (xs)dW ho {x). 


The quantity dW p (x) can be inferred using the QED result: 

dW p ,{x) 


dx 



1 — X + -x z 


L-( 1 -x) 


m n 


\> x = -£-> -±-, L = In 1 + 


ElQ l 

M- 


,P' = P + ,P ,P°, 


(1-9) 


( 1 . 10 ) 


where M , m p , E, and E p are the masses and the energies of the initial proton and the emitted 
p-meson (Laboratory reference frame implied). 

For the probability of hard pion emission we have 


with 


d\V K 

dx 


£|^«(Pi,Pi ~k)Y 


d 3 k 

16(0K 3 : 


I 


& pn (pi,Pi-k)\ 2 =- rr - 

4(P1 % 2 

Hi -2(p x k)} 2 ' 


j—^ Tr (p l -k + M )y 5 (p l +M)j5 = 

( Pl k)=Eco(l-bc),l-b 2 ^’^ + ^, 


( 1 . 11 ) 


( 1 . 12 ) 
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Figure 5: Color online. Quantity dWi,/dx for p- and /r-me son production as a function of the meson energy 
fraction x = Ej t /E, from Eqs. ( 1.10 1.13 ), for E = 15 GeV and for two values of the /i-emission angle: 
0q = 10° for p- meson (black, solid) and for 7T-meson (red, dashed), 0 q = 20° for p- meson (green, dotted) 
and for ;r-meson (blue, dash-dotted). 


with c = cos(&,pi). The angular integration in the region 1 — (6^/ 2) < c < 1 leads to 


dw n ,{x) _ jr_ 

8n 2 


dx 


1 - 


m K 

x 2 E 2 


L + ln 


■ + 


m~ 


d(x) xd{x)M 2 


E-ji 




i ^ ^ "‘ n j/ \ i i m n x 

1>i = Y > T' rf(i) = 1 + «V 


, x = 1—x, n' = 7i + , n ,7r , 


(1.13) 


where g = g ppp = g npp ~ 6 is the strong coupling constant.The quantities dW/,(x)/dx as functions 
of the energy fraction x = E^/E (h = p, n) are drawn in Fig. pj, for £=15 GeV and two different 
values of 6q: 10° and 20°. 

The expressions of the integrated probabilities are: 


wi ' = i‘^ dx= &^ A " L+ ^' 

X? 

AP = I 0 tf)-htf) + ^I 2 tf), BP = -/ 0 (xf)+/i(xf), 
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Figure 6: Renormalized probabilities ■ Pp,n, as function of the incident energy for two values of the 
hadron emission angle. 


A*= l -h(x?y,B* = h(x?), 


(1.14) 


where jcf = £'/'/£' with E\‘ = m h - threshold value of the energy of the detected particle, h = p, n. 
The analytic expressions of the functions 4(z), i = 0,1,2 are presented in Appendix. 

The integrated quantities W/ u h = p,71 can, in general, exceed unity, violating unitarity. To 
restore unitarity, we have to take into account virtual corrections: the vertex for the emission of 
a single pion (charged or neutral) from the proton has to include ’radiative corrections’, which 
account for emission and absorption of any number of virtual pions. For this aim we use the known 
expression for the probability of emission of n "soft" photons in processes of charged particles hard 
interaction, i.e.. the Poisson formula for emission of n soft photons W n = (a"/n\)e~ a (where a is 
the probability of emission of a single soft photon) [14]. 

Note that the probability of emission of ’soft’ neutral pions follows a Poisson distribution, 
which is not the case for the emission of charged pions. Fortunately, in our case, it is sufficient 
to consider the emission of one charged pion at lowest order (the process of one charged pion 
emission) plus any number of real and virtual pion with total charge zero. In such configuration, this 
vertex has the form of the product of the Born probability of emission of a single pion multiplied 
by the Poisson-like factor: 

Pn.p=e- W ^, (1.15) 

which takes into account virtual corrections. 

The renormalized probabilities W^p/Vp < 1 from Eqs. (1.14,1.15) are illustrated in Fig. j|as 
a function of energy, for two different Go angles. 

Keeping in mind the possible processes of emission of n real soft neutral pion escaping the 
detection, the final result can be obtained using the replacement 


k—n 


a{s) -> cr(s) x @ k = P k Y j 


k =o 


k\ 


(1.16) 
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Figure 7: Renormalization factor R K , for the probability of emission of 2 (black, solid line), 3 (red, dashed 
line), 4 (green, dotted line) neutral pions, as a function of the incident energy for 9q = 10°. 


The renormalization factor is illustrated in Fig. [7], for the probability of emission of 2 (black, 
solid line), 3 (red, dashed line), 4 (green, dotted line) pions. 


The quantity P K can be compared to the experimentally measurable phenomena fllC|]: the 
fraction of protons in the final state of proton-proton collisions is approximately one half, P K ~ 0.5. 
The commonly accepted explanation is that charge exchange reactions are responsible for changing 
protons into neutrons. 

Let us consider the antiproton-proton annihilation into two and three pions p + p —>• K + 7i~ or 

p + p —)► ^ + 7r _ 7r°. 

Concerning the production of two charged pions, accompanied by a final state X, we can write: 


da pp ^ p x = 2 


dW p (x) 
dx 


.pp^tX 


(xs) xP p , 


(1.17) 


where the factor of two takes into account two kinematical situations, corresponding to the emission 
along each of the initial particles and P p is the survival factor ( 1.15 ) which takes into account virtual 
radiative corrections. The characteristic peak at x = x max has the same nature as for the QED process 
e + + e~ —>• p + + jU " + y. As explained in Ref. & it is a threshold effect, corresponding to the 
creation of a muon pair, where x max = 1 — 4A7“ /s, A7 i( is the muon mass. 

The cross section ( |1.17[ ) is illustrated in Fig. || for two different values of the laboratory energy 
and of the emitted angle as a function of the p meson energy fraction. 

In case of three pion production,assuming that the process occurs through a n °p° initial state 
emission, we find: 


do(p,p) pp ^ KpX = dWp (x p )dW%{x n )[do{p- p p ,p- p n ) PP ^ X 
do(p- P7l ,p-pp) pp ^ x }P K Pp, 


(1.18) 
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Figure 8: The cross section da{p,p —> p°X) is plotted as function of the p energy fraction for two values of 
the incident energy and of the p emission angle: E = 10 GeV and 0q = 10° (black, solid line), E = 10 GeV 
and Go = 20° (red, dashed line), E = 20 GeV and do = 10° (green, dotted line), E = 20 GeV and 6q = 20° 
(blue, dash-dotted line). 


implying the subsequent decay p° —>• n + n~. 

It is interesting to note that the cross sections for the interaction of high energy neutron (anti¬ 
neutron) beams with a hadronic target can be calculated using the cross sections of proton beam 
interacting with the same target with the emission of the charged meson. We obtain (see Eqs. 

dll): 


_ da pT ^ h+x /dx 
dW + (x) / dx 


(1-19) 


and similarly for the anti-proton beams. 

Experimental data from coincidence experiment with the selection of a hard meson at small 
angles are not available, and we can not test the factorisation prediction. However the values of the 
cross section for ppp°X illustrated in Fig. |] are not in contradiction with the corresponding total 
cross sections (for a compilation, see [16]). 

Inversely, if we use the experimental data for the total cross section of process pp —>• hh p ~ 
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1 mb, we can predict the value of the total cross section of process np —> X 

PnW n (E \, e 0 )a np ^ x (E -E l ) = o pp ~* nX (E), 


( 1 . 20 ) 


with W n (E , 0q) given in Eqs. (1.10,1.16). 


2. Conclusions 

We have extended the QRE method to light meson emission from an (anti)proton beam and 
calculated the probabilities and the relative cross section for multi-pion emission. The considered 
processes can be measured at present and planned hadron facilities. We have also suggested a 
possible application. The collinear light meson emission could also be used to produce secondary 
(anti)neutron beams, at a high energy (anti)proton accelerator. This would constitute an alternative 
to the usual way, when high-energy neutron beams are produced as secondary beams, by break-up 
of deuterons on a hadronic target. 


Our result for the matrix element squared. Eq. ( |1.12| ), agrees with Adler principle: the cross 
section vanishes when the pion momenta vanishes. This result can be inferred as well from the 
reduction formula of current algebra [|n|]. Note that the probabilities to create a n or p-meson by 
a proton, can also be obtained using the infinite momentum reference frame, (Ref. [|l8|], Eq. (52)). 

This allows one to obtain the relation between the cross section dCk(s(l —x)) and dO\{\). 
This relation has the form of conversion on the energy fraction of the emitted particle x = k° { /p ( \ 
of the probability of emission by the initial projectile dWk(x ) with cross section d<7k(s( 1 — jt)). We 
underline that the details of interaction of projectile with the target are irrelevant, assuming that 
they are the same for shifted and unshifted kinematics. The contribution of the final state hard 
hadron emission as well as the interference of the relevant amplitude with the one describing the 
initial state emission vanish in the limit 6q —>• 0. 

Let us discuss the limits of our approach and the accuracy of the present calculation. Con¬ 
cerning the kinematical region where this mechanism is important, we consider charged pion emis¬ 
sion at angles, smaller than an angle do, i.e. 0 < 0o <C 1- Moreover EG » M, which insures 
the presence of large logarithm. In this case, the main contribution is due to the large quantity 
L = \n(E 2 6 2 /M 2 ) 1. In the calculations, we keep the terms containing L and omit the terms 

proportional to 6 2 . The resulting accuracy is given by 1 + 0(1 /L). Taking L = 10, the accuracy of 
the calculation is of the order of 10%. 

We considered the case when the initial and the virtual projectiles are fermions. For instance 


we used the density of probability to find a fermion into the initial fermion (see (E3.E4)) p e e (z) = 
[1 + (1 — z) 2 ]/z- Instead one can consider other relations of the same kind, which hold when the 
initial projectile is a boson (photon, p, meson, gluon). We will not disuss this topic here. Another 
interesting possibility is the conversion of the initial lepton to a photon or a neutral Z boson. In this 


case we must use Pj(z) = zr + (1 — z) 2 [15]. 

The arguments given above have a phenomenological character and are formulated in terms 


of hadrons. A similar idea, at quark level, was introduced in Ref. [19], where the emission of p- 
meson by quark and the p meson production in quark-antiquark annihilation was studied. Special 
attention was paid to polarization phenomena of the created p-meson. 
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We have also suggested a possible application. The collinear light meson emission could also 
be used to produce secondary (anti)neutron beams, at a high energy (anti)proton accelerator. This 
would constitute an alternative to the usual way, when high-energy neutron beams are produced as 
secondary beams, by break-up of deuterons on a hadronic target. 

In frame of the "Gluon Dominance Model", developed by one of us [120]] the ratio of the 
inelastic CE cross section to the total inelastic cross section in pp scattering is estimated as 40%, 
in reasonable agreement with the experimental data [[I]]. 

The collinear light meson emission mechanism in (anti)proton-proton collisions provide a 
possible source of events with rather high multiplicities of (charged and neutral) pion production. 
For such events the emission of hadrons in initial as well as in final states must be taken into 
account. 

We considered here only initial proton emission. In this case, the resonance formation is 
kinematically forbidden since the four momentum of the virtual nucleon is Space-Like. The sim¬ 
plest CE processes pp —>• nil, pn —>• A£ can be in principle measured at PANDA. Other reaction 
mechanisms can contribute to these processes. In the frame of a description in terms of a single 
pseudoscalar meson n + , K exchange, information on the strange meson-baryon constant can be 
extracted. Note that neutrons created by p and n emission form a (anti)proton beam with the mech¬ 
anism considered here (initial state radiation) can not be formed by subsequent decay of a nucleon 
resonance, as A. Nucleon resonances can not be produced in the initial state, because the neutron 
has the time like momentum squared (the propagator is (p — k ) 2 — M 2 < 0) . 

The simplest CE processes pp —>• nh, pn —>• AL~ can be in principle measured at PANDA. 
Other reaction mechanisms can contribute to these processes. In the frame of a description in 
terms of a single pseudoscalar meson K f , K 1 exchange, information on the strange meson-baryon 
constant can be extracted. 

3. Acknowledgements 

One of us EAK is grateful to RFBR grant 11-02-00112 for support. We are grateful to S. 
Barkanova, A. Alexeev, and V. Zykunov, for interest to this problem and to V.A. Nikitin for useful 
discussions. 


4. Appendix 

The analytic expressions for calculating the integrals 


are: 


•w=/H 


W) = 

Hz) 
h ( z) 
Hz) 


i-(- 

x 


1 1 +r 

— In-r; 

2 1 - r 

r + zarcsin(z); 


-r -In - 

2 4 l—r 


with r = \/l — z 2 - 
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5. Instead of Conclusions: A memory of Professor Eduard Alekseevich KURAEV 

Prof. Eduard Alekseevich KURAEV passed away, while working on this subject. 

It is an invaluable loss for our community, in particular for the colleagues who had the chance 
to meet him and to share his encyclopedical knowledge, his enthusiasm, and his original and 
surprising ideas. 

Always concentrated on solving difficult problems. He published hundreds of scientific papers, 
and collected thousands of citations. 

He had an astonishing energy, extraordinary ability and resistance for effort and hard work., 
a passion for sharing ideas, discoveries, and results in seminars and in teaching and educating 
young students and researchers. He was always ready to listen and discuss with the colleagues 
knocking at his door, coming from near or far Institutes. 

He had a deep understanding and a curiosity in the observation of Nature, of the visible and 
invisible features of Nature, and, at the same time, a strong capability in abstraction and a passion 
for mathematical problems. 

We loose with Prof. Eduard Alekseevich Kuraev, an engaging and unique personality, and an 
exceptional physicist. 

References 

[1] V. S. Murzin and L, I. Sarycheva. Multiparticle processes at high energies. Atomizdat, M. (in Russian) 
1974, p. 253, and references therein; S. Baskovic et al. Yad. Phys., 27, 1225 (1978). 

[2] V.V. Avdeichikov et al. Instrum. Exp. Tech. 56 (2013) 9. 

[3] M. A. Braun, Y. .A. Evlashev and V. VVechernin, Yad. Fiz. 27 (1978) 1329. 

[4] E. A. Kuraev, S. Bakmaev and E. A. Kokoulina, Nucl. Phys. B 851 (2011) 551. 

[5] E. N. Antonov, L. N. Lipatov, E. A. Kuraev and I. O. Cherednikov, Nucl. Phys. B 721 (2005) 111. 

[6] L. N. Lipatov, Nucl. Phys. B 365 (1991) 614. 

[7] S. Gasiorowicz, Elementary particle physics, J. Wiley & Sons Inc., New York (1966). 

[8] V. Baier, V. Fadin and V. Khoze, Nucl. Phys. B65, 381 (1973). 

[9] A. B. Arbuzov, V. V. Bytev, E. A. Kuraev, E. Tomasi-Gustafsson and Y. .M. Bystritskiy, Phys. Part. 
Nucl. 41, 593 (2010); Phys. Part. Nucl. 42, 1 (2011). 

[10] V. A. Nikitin, private communication. 

[11] M. Caffo, R. Gatto and E. Remiddi, Nucl. Phys. B 252 (1985) 378. 

[12] Physics Performance Report for PANDA: Strong Interaction Studies with Antiprotons, The PANDA 
Collaboration, arXiv:0903.3905 [hep-ex]; http://www.gsi.de/PANDA; 

[13] http://www.gsi.de/FAIR; 

[14] A. I. Akhiezer and V. B. Berestetskij, Quantum Electrodynamics, Moscow Science, 1981. 

[15] V.N. Baier, V.S. Fadin and V.M. Katkov, Emission of Relativistic Electrons, Atomizdat, Moscow, 

1973 (in Russian). 

[16] A. Dbeyssi and E. Tomasi-Gustafsson, Prob. Atomic Sci. Technol. 2012N1, 84 (2012). 


13 




Light meson emission 


E. Tomasi-Gustafsson 


[17] I.I Sakurai, Chicago Lectures in Physics, The University of Chicago Press Chicago and London 
(1967) Chapter 5. 

[18] G. Altarelli and G. Parisi, Nucl. Phys. B 126, 298 (1977). 

[19] A. V. Efremov and O. V. Teryaev, JINR Preprint P2-82-832 (in Russian). 

[20] E.S. Kokoulina, Acta Physica Polonica. B35, 295 (2004); M.V. Volkov and E.S Kokoulina, Phys.Part. 
Nucl. Lett. 1, 16 (2004). 


14 




